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We show that LHC experiments might well be able to determine all the parameters 
Q^' required for a prediction of the present density of thermal LSP relics from the Big Bang 

era. If the LSP is an almost pure bino we usually only need to determine its mass and the 
D , mass of the SU{2) singlet sleptons. This information can be obtained by reconstructing 

\ the cascade ql X2Q ~^ ^R^Q Xi^~^^~Q- The only requirement is that m^^ < m^o, 

which is true for most of the cosmologically interesting parameter space. If the LSP has 
^ . a significant higgsino component, its predicted thermal relic density is smaller than for 

^ I an equal-mass bino. We show that in this case squark decays also produce significant 

numbers of Xa aiid xt'- Reconstructing the corresponding decay cascades then allows to 
determine the higgsino component of the LSP. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) |I| is one of the most promising 
extensions of the Standard Model. It offers a natural solution of the hierarchy problem 
amazing gauge coupling unification 0], and Dark Matter candidates 0]. If Nature chooses 
low energy supersymmetry (SUSY), sparticles will be found for sure, as they will be copiously 
produced at future colliders such as the Large Hadron Collider (LHC) at CERN or TeV scale 
e~^e~ linear colliders (LC). The LHC would be a great discovery machine if SUSY breaking 
masses lie below a few TeV [0. On the other hand, there are several on-going and future 
projects searching for LSP Dark Matter. One of them even claims a positive signal @], although 
the current situation is rather contradictory . In any case, it seems very plausible that both 
SUSY collider signals and LSP Dark Matter in the Universe will be found in future. 

Recently interesting possibilities have been pointed out where non-thermal production of 
Dark Matter is significant P, Generally the known bound from the thermal LSP density 
may easily be evaded by assuming a low post-inflationary reheating temperature of the Uni- 



verse, without endangering the standard successes of Big Bang cosmology |]T0| . If the reheating 
temperature is below the neutralino decoupling temperature, the relation between neutralino 
pair annihilation rates and the mass density of the Universe disappears. The mass may then 
be determined by other parameters, such as the Q ball formation rate and decay time 0, or 
the moduli masses and their decay rates to LSPs 0. 

While these non-thermal mechanisms open exciting new possibilities, direct experimental 
or observational tests of them might be difficult, since they all have to occur before Big Bang 
nucleosynthesis (BBN).[]It is therefore interesting to determine 

1. the actual LSP relic density, both "locally" (in the solar system) and averaged over the 
Universe; and 

2. the predicted thermal LSP relic density, 

as precisely as possible. These quantities are closely related to the mass and interactions of 
the LSP. A positive difference between the actual and predicted LSP density would indicate 
the existence of non-thermal relics, whereas a negative difference would hint at large entropy 
production below the LSP freeze-out temperature (e.g. due to a low reheating temperature). 

The matter density in the Universe divided by the critical density, ^matter, is claimed to 
be tightly constrained already; ^matter = 0.35 ± 0.07 |jl2|. On the other hand, the thermal 
relic density of the Universe ^^oh^ {h = 0.65 ± 0.05) has been calculated through the mass 
and interaction of the LSP, which is likely to be the lightest neutralino Xv the absence 
of direct production of sparticles we have to rely on experimental lower bounds on sparticle 
masses as well as naturalness arguments to conclude that the predicted ^j^oh? lies somewhere 
between 10~^ and 10^; clearly this is not a very useful prediction, although it is encouraging 
that this wide range at least includes the correct value. The purpose of this paper is to discuss 
how future LHC experiments can contribute to the determination of the MSSM parameters 
that are needed to predict the thermal LSP relic density and the LSP-nucleon scattering cross 
section. Our goal is thus somewhat different from that of ref . [|l^ , where it was simply assumed 
that all relevant parameters had somehow been determined by various experiments, with given 
errors; the main emphasis was on estimating the resulting uncertainties in the predictions of 
the thermal LSP relic density and the LSP-nucleon scattering cross section. In contrast, we 
discuss in some detail how these parameters can be determined, and with what errors. 



-"^For the Q ball case, measurements of the cosmic microwave background by the MAP and PLANCK 
satellites might find isocurvature fluctuations due to the Affleck-Dine Field 
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The determination of mass parameters has been discussed in detail in the minimal super- 
gravity (mSUGRA) model, where one assumes universality of scalar masses and of gaugino 
masses at the scale Mx of Grand Unification In Sec. 2, we point out that X2 ^ is 
open for most parameters giving a reasonable LSP density, making the determination of m^^o , 
m^o and m^^ possible at the LHC. We demonstrate that the mass density is determined by 
the LSP and slepton masses, if the LSP is mostly a bino as expected in mSUGRA. In this 
case fi-j^o can be predicted to about 10 to 20% accuracy. 

In Sees. 3 and 4 we discuss a non-mSUGRA scenario. In Sec. 3 we relax the assumption 
of universal scalar masses for Higgs bosons. It is then easy to find cases with comparable 
higgsino and gaugino masses, ~ M, while keeping all squared scalar masses positive at Mx- 
The LSP then has a significant higgsino component, so that its density cannot be predicted 
by only studying X2 ~^ ^ Xi decays. The situation is further complicated if we also relax 
the assumption of universal gaugino masses, since the neutralino mass matrix then depends 
on three independent, unknown mass parameters. We point out that the cascade decay xt ~^ 
ul — >■ Xi can then often be identified, providing clear evidence that /i ~ M. In Sec. 4 we 
present a detailed case study with /i ~ M2 to confirm the potential of LHC experiments to 
analyze xt cascade decays; this allows a complete determination of the neutralino mass matrix 
(in the absence of CP-violating phases). Sec. 5 is devoted to discussions. 



2 Q~o in mSUGRA 

Al 

In the minimal supergravity model one assumes universal soft breaking parameters at the GUT 
scale: a universal scalar mass m, universal gaugino mass M, universal trilinear coupling A, 
and Higgs mass parameter B. The renormalization group evolution of soft breaking squared 
Higgs masses then leads to consistent breaking of the electroweak symmetry, provided the 
higgsino mass parameter /i can be tuned independently. In this paper, we chose the weak 
scale input parameters mh{mh) = 4.2 GeV, mt{mt) = 165 GeV, and tan/3. We minimize the 
tree level potential at renormalization scale Q = ^ymfrn^, which reproduces the correct value 
of fi obtained by minimizing the full 1-loop effective potential. We include loop corrections 
to the masses of neutral Higgs bosons, including leading two-loop corrections ||15||. 



The mass density fi^?^^ '^^ IjSP is calculated from the pair annihilation cross section 



by using the expressions [|I 



1.07 X lOyCeVxp 
y/g^Mp{a + 3b/ xp) 
aixiXi ^ ^^^)v = a + bv^ 



Al 



xp{= m^o/Tp) = log 



0.764Mp(a + 6b/xp)c{2 + c)m^o 



except near regions of parameter space where special care is needed. Here Mp = 1.22 ■ 10^^ 
GeV is the Planck mass, and a and b are the first two coefficients in the Taylor expansion of the 
pair annihilation cross section of the LSP with respect to the relative velocity v of the LSP pair 
in its center of mass frame, g* is the effective number of relativistic degrees of freedom at LSP 
freeze-out temperature Tp. The expansion in v breaks down around s— channel poles; here 
the thermal average is calculated numerically, using the formalism given by Griest and Seckel 



1^. We also take into account sub-threshold annihilation into hh and W^W final states. 
When the LSP is higgsino-like, coannihilation of x? with xt ^"^^ X2 ^'^c important [18, 19 . 
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The annihilation modes XiXt ~^ ff X1X2 ~^ ff approximated by s— channel W and 
Z exchange, respectively; XiXi ~^ ^1 is also included. All other higgsino coannihilation 
modes are treated assuming SU (2) invariance. We also include one loop radiative corrections 
to the mass splitting of higgsino-like states pO| . We do not include Xi^ coannihilation, since 
we do not study cases with m^o ~ m^. 

MSUGRA predicts a bino-like LSP Xi and wino-like Xi and X2 f*^^ moderate values of m 



and M (below ~ 500 GeV). This is a rather model independent result [^. Large positive 
corrections to squark masses from gaugino loops, together with the large top Yukawa coupling, 
drive the squared soft breaking Higgs mass m2 negative at the weak scale. On the other hand, 
correct symmetry breaking requires ml + fi'^ > — m|/2. One has to make /i large to obtain 
the correct electroweak symmetry breaking scale, if scalar masses and gaugino masses are of 
the same order. 

If slepton masses are moderate, the LSP is bino-like, and one is sufficiently far away from 
s— channel poles (2m^o 7^ mz, "n^Higgs), the mass density is essentially determined by t— channel 
cr exchange This is because 



1. A pure bino couples only to fermion and sfermion, or Higgs and higgsino. Higgsino 
exchange is suppressed for /i^ ^ M^. 

2. mg^ < rriej^ — mc, <^ rriq in mSUGRA, therefore cr exchange is least suppressed by large 
masses in the propagator. 

3. The hypercharges of sleptons satisfy the relation Y^^ = 2Ye^, therefore a{eR)v = 16a{eL or h'L)v 
when sfermion masses are equal. 

This can be seen in Fig. 1 a)-d), where fl^o (a, b) and b = 10® GeV^ x Q^ah'^aj^ (c,d) 
are plotted. Here o"^ is the scaled bino pair annihilation cross section in the limit where 



rrier, <^ mg. , mg 19 



a- - ^° 



'B 



m? + m-o / (rni + m-n)^ 



(2) 



We find that the mass density increases with increasing M oc m^o and m; m is essentially 
proportional to m^^ for m > M. Dotted lines are for constant Xi mass. The fact that it 
basically only depends on M indicates that the LSP is indeed bino-like. The mass density 
becomes very small for m^o = mz/2, because LSP pair annihilation through Z exchange is 
enhanced. In Fig. 1 c) and Id) we show contours of constant b. Although VL^oh? changes by 
more than a factor of 4, the change of b is very small over the wide range of parameter region 
with M > 160 GeV, again confirming the bino-like nature of the LSP for the mSUGRA case. 

The tan/3 dependence is also very weak, as can be seen in Fig. 2 a) and b). This again 
shows that the LSP is bino dominant, and bino-higgsino mixing, which is controlled by the off- 
diagonal elements of the neutralino mass matrix, has negligible effect on the LSP relic density 
for the parameters given in the plot. For sufficiently heavy ^x? simply determined by 

rrieji ^^"i '^x? ^^^^ 6 ~ 1. 

We now show that analyses of sparticle production at the LHC would lead to tight con- 
straints on the predicted thermal relic density fl^oh"^. Recently, quite a few studies of precision 
measurements of sparticle masses at the LHC have been performed. When the cascade decay 
g ^ X2 ^ ^-R ^ Xi is open, a clean SUSY signal is // + jets+ missing Et- It was shown 
that rriq, m^o , and m^o can be reconstructed from the upper end points of the mju and mji 
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Figure 1: Contours of constant ^j^ab? (Fig 1 a,b) and h factor (Fig 1 c,d) in the (m,M) plane for tan/3 = 10 
(a,c) and 4 (b,d). We take /x > 0. Contours of constant ?n^o, as well as contours where m^a = mgR and 
m^o = rrie^ are also shown. 
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distributions, mjlf^ and m™^^; the edge of the mu distribution, m™'^^; and the lower end point 
of the rrijii distribution with ma > m^^^/\/2, mjl]^. Here j refers to one of the two hardest 
jets in the event. In most cases it is chosen such that it has the smaller jU invariant mass; this 
is meant to select the jet from the primary q X2I decay. However, m™™ is reconstructed 
by taking the jet which gives the larger jll invariant mass, in order to avoid contamination. 



Those end points are given by the analytical formulae [14 



m 



m 



m 



'jll 



(m|^-m|o)(m|o-m|o) 



1/2 



Max 



V 



m-o) (m-o — m? 1 



1/2 



(m? — m?o)(m? — m?o 



2 

mi 



1/2^ 



(m|„ -m|o) 



4m|omi^ 

2 2 
— m-om 



-m|om|o + 3m|om|om?^ - mjom^^ - m|om|^ - m|om|om?^ 



^ 4 2 



y(m|o+m|^)(m|o+m?^; 



2 + 2m-om 



efl^ X^ 



■ Qminml + m^ 

X2 ^fl' 



(3) 



In addition to those quantities, one can measure the end point m^™ of the distribution of 
the smaller of the two m^; values. It can be expressed as 



m 



\ 



\ 



if 2m' 



2mi^ 



m?o 

X? 



(m|o + m|o) < 

^ Xi X2' 



(m? —m-o){m-o—m^ ) 

it 2m. - (m-o + m-o) > 

^Jr<. ^ Al A.2 



m?o 

xs 



(4) 



Because there are only four masses involved, the last end point is redundant, but might be 
useful to cross check the decay kinematics.0 

For the example studied in |jT^, the so-called "point 5", m = 100 GeV and M = 300 GeV, 



which results in m^o 

X2 



233 GeV, m 



157.2 GeV and m^o 

Xi 



121.5 GeV. The errors on 



and m^o are strongly correlated and are found to be 12% for m^o and 9% for mg^. Within 
the framework of mSUGRA the measured LSP mass excludes the possibility that s— channel 
poles are important for the LSP pair annihilation cross section (see below). We find that the 
corresponding error on cr^, and hence on the prediction for fl^oh'^, for this parameter point is 
20%. If the error (which is dominated by systematics associated with uncertainties of signal 
distributions) is reduced by a detailed study of various signal distributions, the error on o"^ 
may go down below 10%. 

^We are assuming that squarks are basically degenerate. Note that essentially only left-handed squarks 
will contribute here, since SU{2) singlet squarks very rarely decay into X2- Barring "accidental" cancellations, 
bounds on flavor changing neutral current processes imply that squarks with equal gauge quantum numbers 
must be close in mass. The mass splitting between ul and c^l squarks is limited by SU{2) invariance. With 
the possible exception of third generation squarks the assumed degeneracy therefore holds almost model- 
independently. If required, contributions from third generation squarks can be filtered out by anti-tagging 
6— jets. 
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Fig. 1 also contains contours where m^o = mg^ and m^o = mg^. To the left of these 

contours X2 decays into / are accessible, giving a substantial constraint to the kinematics of 
the events. First, note that X2 ~^ decays are open for most of the cosmologically acceptable 
region with M > 200 GeV. There is even a substantial region of parameter space where both 
X2 ~^ and X2 ~^ are open. In our argument above, we assumed that all observed edges 
and end points kinematic distributions come from Ql ^ X2 ^ ^R rather than X2 ^ ^^l- When 
the latter decay mode is open, the branching ratio dominates over eji, since X2 — ^3 and 
is an SU{2) singlet. However if one assumes that squared scalar masses are positive at the 
GUT scale, cannot be too much lighter than x% so there is some kinematical suppression. 
In mSUGRA the relevant masses are expressed as: 

Ml = 0.4M, M2 = 0.8M 



+ O.SM^ ~ ~ ^^^^ 
ml^ = + 0.15 -sm^Owml cos 2p, (5) 

and m^o ~ Mi and m^o ~ M2 if M ^ /i. In it was shown that the two edges can be 
observed separately even if X2 ei is not strongly phase space suppressed. It might also be 
possible to find evidence for light left-handed sleptons by looking into the relative strengths 
of different SUSY signals. If X2 ~* z/^ is open, Xi ^ (^l or is also open and dominant, 
yielding relatively large and signals compared to the signal. 

If we assume bottom Yukawa corrections are negligible and squared scalar masses are 
positive at the GUT scale, the pseudoscalar Higgs mass is bounded from below as m\ > 



/i Under the bino dominant assumption, and for moderate value of tan/?, neutralino 



annihilation through s— channel poles can thus not be important. On the other hand, for 
large tan/? the pseudoscalar Higgs boson can be light enough to achieve m^o ~ mA/2 p2|, p!9[| . 
However, in mSUGRA large tan/? also implies a rather light fi, which greatly depletes the 



/"•"/" signal ||2^; the observation of strong multi-lepton signals would thus already indicate 
that tan/5 is not very large. Note also that for large tan/? direct production of the heavy 
neutral Higgs bosons from gluon fusion and/or in association with bb pairs allows to detect or 
exclude these Higgs bosons at the LHC for up to several hundred GeV . We will come 
back to the importance of determining later in Sec. 3. 

In Fig. 1 and 2, we only looked into the parameter space with moderate ("natural") values 



of m and M. If m ^ M, solutions with /i ~ M may be obtained |]T9|, |2^ , and the assumption 
X? — -B is not valid any more. In such a case the decay X2 ~^ ^R '^'^^ open. However 
there is still a chance that wino-like charginos and neutralinos (^4 and xt) are produced in 
cascade decays, and yield additional kinematic constraints besides the end point measurement 
of the mil distribution from x% — > xV-l decay; see Sec. 3. Large m may also be allowed when 
m^o ~ vtizjl, in which case X2 ~^ ^^R need not be open to make f2^o small. Even then the mu 
end point determines m^^o — rrij^o [|, and the mu distribution of three body X2 decays is 



sensitive to very large mg^ Another twist appears when X2 ~^ Xi^ is open and dominates 
X2 decays. The small leptonic branching ratio of the Z boson might then make it difficult to 
study neutralino masses, and there is no sensitivity to slepton masses. Note, however, that 
in the bino-dominant region Z decays cannot compete with X2 ~^ ^R decays unless the 

latter are strongly phase space suppressed. 

We already briefly alluded to the case where ^ m^^ due to renormalization group 
effects and f mixing. The lighter f can be substantially lighter than the other sleptons if 
/itan/5 is large |2^. In this case pair annihilation through t channel f exchange can even 



dominate other sparticle exchange contributions [1^, because ri could be lighter than the 
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other sparticles, and the mixing induces an S"— wave amphtude. In p8| the possibihty to detect 
and study f at the LHC is discussed. The end point of the jrjr invariant mass distribution, 
where jr denotes a r— jet, is not as well determined as that of the mu distribution, but it 
has been estimated that a ~ 5% measurement should be possible. Even if the jVjV end point 
indicates m^^ < rrie^j, the constraint on rrig, m^o and m^o from // events originating from 
X2 ~^ decays can perhaps be used to reduce the mf error, in which case the combined error 
should not increase too much. However, if tan/3 is very large, it becomes easier to have an 
acceptable LSP relic density even if mg^ > frtxl' ^^^^ ^^^^ ^^^^ need a linear collider 
to perform precision measurements of the nature of fi |2^, where 0"^^, the end point of 
the r jet energy distribution, and a measurement of the r polarization would do a good job 
in determining the parameters needed to predict the thermal LSP relic density. 



3 in non-mSUGRA scenarios and collider signals 

In the previous section, we have shown that the mSUGRA assumption predicts a bino- 
dominant LSP. We also found that measurements at LHC experiments are sufficient for a 
prediction of VL^^h'^, if the cascade decay X2 ^ ^ Xi is open and LSP bino dominance 
is assumed. Now the question is if LHC experiments can be used to check the assumption 
that the LSP is mostly a bino. After all, it is possible that is smaller than or of the order 
of the gaugino masses, and that the mSUGRA relation Mi ~ M2/2 is broken. In this and 
the following Section, we discuss a scenario where the inequality Mi < M2 is kept, while 
/i is substantially smaller than the mSUGRA prediction. In such a case Z exchange effects 
and/or LSP annihilation into W pairs are expected to be more important than in the strict 
mSUGRA scenario studied in the previous Section, and one needs more information to predict 
the thermal contribution to VL^oh?. 

The relative size between /i and M is controlled by Higgs sector mass parameters. The 
MSSM Higgs potential can be written as 

V ={ml + fi^)HlHi + {ml + /i^)i/li^2 + {B^iHi ■ H2 + h.c.) + (4th order terms). (6) 

Here mi and m2 are soft breaking Higgs masses. In the previous Section we took mi = m2 = m 
at the GUT scale, which gave /i^ ^ M^ unless th? 3> M^. In general, ~ M may be achieved 
by allowing non- universal soft breaking Higgs masses, mi, m2 7^ m. For simplicity we will 
keep mi = m2 = m^ at the GUT scale; we briefly comment on the effect of relaxing this 
assumption below. 

In Fig 3a), we plot vs. (m/i/m)^. By increasing rrih, is reduced gradually so that 
m\ + at the GUT scale is roughly constant. Note that the negative radiative correction to 
Higgs mass m\ is dominated by the gaugino mass through the stop and sbottom masses for 
this choice of parameter^. Therefore the value of m\ giving |/i| ~ M2 is almost independent 
of m, rrih ~ 330 GeV. Generally ~ M can be achieved if mh\GUT ~ the precise value 
is determined by the top Yukawa coupling. In Fig. 3b) and 3c), we also plot f2^o and the h 
factor. These quantities vary substantially once falls below M2. 

In Fig. 4 we show contours of constant VL^o and constant h factor in the {m\/Tn?, M) plane 
for flxed m = 100 GeV. We flrst note that fl^o decreases as ruh is increased. This is due to the 
reduction of for larger m/j, see Fig. 3. This increases the higgsino components of Xi and 

■^The fact that in mSUGRA at the weak scale is ahiiost independent of is closely related to the 
"focus point" behavior studied in ref.p^. 



7 




Figure 3: a.) ijh) ^l^o and c) b as function of {mh/mf. We fix M = 250 GeV, A = and tan/3 = 10. 
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also reduces its mass. Not only f2^o but also the b factor decreases, therefore pair annihilation 
is no longer dominated by sfermion exchanges. Especially in Fig. 4d) we find very strong 
effects from LSP annihilation into W pairs. The same effect also can be found in Fig. 3b and 
c), where the rise of b and Q^o corresponds to the closure of the WW mode. Note that no 
consistent solution with electroweak symmetry breaking exists below and to the right of the 
dashed line.0 

As mentioned earlier we assume the two soft breaking Higgs masses to be the same at the 
GUT scale. However, once tan^/3 ^ 1, the higgsino mass is essentially only sensitive to 
the value of rn^. One can therefore increase or reduce the pseudoscalar Higgs mass by varying 
ml at the GUT scale without affecting /i at the weak scale significantly. Nevertheless, as long 



as ml{Mx) > and tan/5 is not very large, rriA = y^i + ml + 2fi'^\^eak remains well above 
2m^o. However, in principle there is nothing wrong with having m^(Mx) < 0, as long as the 
"boundedness condition" ml + > 2|i?/i| remains satisfied at all scales. If mA — 2m-^o, the 
thermal LSP relic density is very small. Strictly speaking the constraints on VL^oh? that we 
will derive below are therefore merely upper bounds as long as we cannot prove experimentally 
that mA is well above 2m^o. 

The reduction of would alter SUSY signals at colliders significantly. When < M, x\ 
and xt production from the decay of SU{2) doublet squarks becomes important as they have 
substantial wino component. This leaves an imprint on the kinematics of di-lepton events, 
which gives us access to additional MSSM parameters, especially when the decay channels of 
neutralinos and charginos into real sleptons are open. This increases the statistics of clean // 
+jets+missing Pt events, since the channels 

~+ -{*) ~-|_ 

Xl - et)(eW)-x^(x?) (7) 

should be seen in addition to the conventional X2 Xi^^ signal.0 
As a result, 

1. The mil edge, and the other end points of the jl and jll invariant mass distributions of 
X2 ^ Xi be measured. By identifying r— leptons from xf decay, one can 
confirm experimentally that the cascade decay originates from xt- This gives lower and 
upper bounds on which in turn constrain the size of the higgsino component of the 
LSP. An explicit example will be analyzed in Sec. 4. 

2. xt '^^y decay into ei directly followed by cl decay into X2 or Xi- O^i the other hand the 
decay X2 ^ is usually forbidden or kinematically suppressed. Therefore X2 decays 
and X4 decays give information on different slepton masses. 

Note that there are substantial constraints on -ino masses and slepton masses from SU (2) x 
f/(l) gauge invariance. The six chargino and neutralino masses are determined (up to radiative 

^Right on this line electroweak symmetry breaking requires ^ = 0. Searches for neutralino and chargino 
production at LEP therefore exclude the region just above the dashed line. However, this experimentally 
excluded region is very narrow, since varies very rapidly near the maximal allowed value of mf^ as shown 
in Fig. 3. 

^ When Ml and M2 have the same sign and is not too small, one of the neutralinos is very higgsino-like 
and would not be produced from the first and second generation squark decays. Here we implicitly assume 
Ml < n ^ M2, in which case the higgsino-like state is Xa- 
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corrections by the values of the four parameters Mi, M2, fi and tan/5, while mg^ and mj> 
are related by 

m| — m^^ = m| cos^ 6'i4/ cos 2/3 (8) 

Therefore the measured edges and end points originating from several decay chains can over- 
constrain the relevant MSSM parameters. 




Figure 5: Squark decay branching ratios for a) m = 90 GeV and b) to = 120 GeV, for M = 250 GeV, A = 
and tan /5 = 10. We average ul and decay branching ratios. 



In Fig. 5, we show various ql decay branching ratios, defined as an average of ul and 
di branching ratios. As ruh increases we find substantial branching ratios into the heavier 
neutralino and chargino, once becomes comparable to M2. The sources of the rise of the 
of X2 signal are: 

1. The increase of squark branching ratios into xt Xa with increasing wino component. 
In the limit |/i| <^ M2 the branching ratios satisfy (assuming Mi < M2): 

Br{qL ^ xt) ■■ Br{qL ^ x1) ■ Br{qL ^ x^) = 2 : 1 : 0. (9) 

Throughout Fig. 5, Xs is a nearly pure higgsino and does not have a substantial branching 
ratio. See Table 2 in Sec. 4 for an explicit example. 

2. The growth of the xt decay branching ratio into I'lh again due to the increase of its 
wino component. Note that xt ~^ is kinematically favored compared to xt ~* (^l- 

3. There is little or no phase space for decays into xtl if M2 < On the other hand, 
this mode may open up when xf becomes higgsino-like. Indeed for Fig. 5a), in the 
region of small the decay ul Xi is suppressed due to the near mass degeneracy 
between ul and xt- 

4. Finally, "conventional" modes such a.s q —>■ xt and q ^ X2 "^iii be suppressed due to the 
reduced gaugino components of these light states. Moreover, the decay X2 ~^ be 
suppressed or closed kinematically even for small values of m, if < M2. In the limit 

< Ml, m^o — m-j^o is small, and the corresponding uiu distribution might be too soft 
to be accessible experimentally. 
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These observations tell us that one should look for X2 and X4 signals in addition to conven- 
tional X2 ~^ Xi^^^~ decays. Discriminating experimentally between scenarios with > M, 
where these new signals are small, and \fi\ < M2, where they are expected to be significant, 
would be important to predict the mass density of the Universe. In the next Section we 
illustrate how these new signals could be analyzed at the LHC 



4 Analyzing the MSSM with \fi\ - M at the LHC 
4.1 A Monte Carlo Study 



m 


90.0 


tan (3 


10 


nil 


360 


1712 


360 


M 


250 


fi 


199.85 


Ml 


103.9 


M2 


208.75 




139.3 


mg^ 


206.09 




556.07 




561.67 




190.28 




93.18 




534.36 




533.21 




155.13 


A3 


208.74 




148.44 




272.52 




273.8 




188.67 




374.43 




563.81 




132.56 




206.06 




498.27 




531.2 


rrih 


112.59 


mp 


436.98 


niH 


437.63 




624.36 



Table 1: Mass parameters and relevant sparticle masses in GeV for the point studied in this paper. ISAJET 
fs^ was used to generate this spectrum. 

We now study leptonic SUSY signals at the LHC for a case where xt production from qi 



decays is sufficiently common to be detectable. We used ISAJET 7.42 to generate signal 
events, while ATLFAST 2.21 was used to simulate the detector response. For this analysis 
we choose the MSSM parameter point shown in Table 1. Here we took a moderate value for 
M, leading to a large sample of signal events. 

The value of the GUT scale Higgs mass is chosen such that the B component of Xi ~ 
0.9, so that effects from its other components on the predicted LSP relic density start to 
be significant; {N^,N^, Njj^,NfjJ = (0.91,-0.15, 0.19,-0.35). In ISAJET this requires 
rrih/rn = 4. Our one-loop RG analysis described in Sec. 3 reproduces this point for rrih/m ~ 
3.74, see Fig. 3a. See also Figs. 3b) and 3c) for the corresponding values of fl^o and the b 
factor. Reducing |/i| even further (by increasing m/j) would lead to even larger differences to 
well-known mSUGRA scenarios, making it easier to measure all relevant parameters. 

In pp collisions one mostly produces squarks and gluinos. They decay further into neutrali- 
nos and charginos. In our case gluinos do not decay exclusively into third generation squarks; 
the branching ratio into first and second generation left handed squarks is about 20% (11.3% 
into ul and 9.6% into di)- These squarks then often decay into xt and xi- The branching 
ratios relevant for the following discussions are summarized in Table 2. The tiny branching 
ratio into Xg is due to the fact that it is mostly higgsino. In Table 3 we show the dominant 
cascade decay processes which produce opposite sign same fiavor lepton pairs in the final state. 
In this Table we also list the corresponding end points of the kinematic distributions discussed 
in Sec. 2, see eqs.(|) and (^). 

We now show several SUSY event distributions, after applying the following cuts |^ to 
reduce the SM background to a negligible level: 
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dl X4 


15.4 


ul xt 


21.4 


dL xt 


36.4 
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23.6 


Xl eL 


5.5 


X2 
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xl e/j 
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Xa n 


89.1 


i'L Xl 


40.1 


CL X2 
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en Xl 


100 


gl Xi 
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Table 2: Relevant branching ratios in % for our sample parameters. 



mode 


lyy-f VYlSiX 

''''jl 


''''jl 


'''jll 


'''jll 


'''ll 


Dl) QL ^X'2^ en ^ X? 


400.3 


237.0 


430.6 


164.2 


50.6 


D2) qL ^xt^i>L^xt 


351.0 


260.0 


411.1 


211.1 


121.9 


D3) QL ^ x2 ^ ^ x^ 


322.7 


269.5 


403.3 


207.5 


117.9 


D4) qL^xl^eL^ Xi 


437.0 


321.4 


460.6 


246.6 


159.6 


D5) Ql^Xa-' en ^ Xi 


421.5 


292.2 


460.6 


270.3 


174.3 



Table 3: End points or edges of invariant mass distributions (in GeV) for different decay processes. 

• 4 jets with Pt,i > 100 GeV and Pt,2,3,4 > 50 GeV. 

• Meff = Pta + Pt,2 + Pt,3 + Pta + > 400 GeV. 

• > Max(100GeV,0.2Mefr). 

• Two isolated leptons with P^ > 10 GeV and |?7| < 2.5. Isolation is defined as having less 
than 10 GeV energy deposited in a cone with AR = 0.2 around the lepton direction. 

In the following plots we reduce SUSY backgrounds by subtracting event samples with 
different fiavor, opposite sign dileptons (e+/i~ and /i~e+ ) from the sum of the e~^e~ and 
event samples. To do this consistently we require two and only two isolated leptons in the 
final state.0 We generated events corresponding to an integrated luminosity of 200 fb~^, but 
the figures are normalized to 100 fb~^. 

In Fig. 6 we show the di~lepton invariant mass distribution for our representative point. 
After the subtraction of e/i events, we see a distribution with at least four edges.[] They are 
consistent with those found in Table 3. Note that a rather weak edge from decay D3) should 
appear very close to the one from D2) in both the M ^ \fi\ and ^ M limit; m-+ ~ ^^x"' 

^However, as we have discussed in previous Sections, the rates of 4 and 3 lepton events compared to 2 
lepton events must contain important information about MSSM parameters. 

^Note that at mu ^ 55 GeV, this subtraction reduces the number of events by a factor of 0.35. The 
fluctuation of the resulting distribution is therefore higher than what is expected from the number of events 
in this distribution. 
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distribution for the parameter point Usted in Table 1. 



m-fO ~ m-+ 

Al 



and mc, ~ mg^ hold in a wide region of parameter space. The two edges must 
be separated out by fitting the smeared mu distribution. Note that since the kinematics of 
two decay chains is expected to be similar, the systematic errors associated with the fitting 
should be small. It seems that at least the first four mu edges can be used for the fit of MSSM 
parameters, while it is not clear if the last one is detectable statistically. 

We then follow the analysis of [Q, by taking the jets with the first and the second largest 
Pt and considering their mju distributions. We label ji and j2 so that mj^^u < rrij^ii. We 
then find that most events have mj^^n below ~ 400 GeV. The mju distribution will contain 
events from the different decay chains listed in Table 3, but they can easily be separated out 
by requiring ma to lie between certain values. For example, if we require that mu < 55 GeV 
(Fig. 7a) and 55 < mu < 125 GeV (Fig. 7b), the distributions should dominantly contain 
events from decay chains Dl) and D2), respectively. In Fig. 7a) the mju end points are indeed 
consistent with the values of end points listed in Table 3. The distribution in Fig. 7b) is 
somewhat smeared out near the end point, due to contamination from Xa decays. 

In the next step we select events where mj^u < 500 GeV < mj^u, the resulting mju distri- 
butions are shown in Figs. 7c) and 7d). These additional cuts have been applied in ref.[0 
because they reduce the probability to select the "wrong" jet, which does not come from 
primary decays. The mj^u distribution is then substantially harder, better reflecting the 
distribution of the "correct" jet. Especially for events with 55 GeV< mu < 125 GeV, the mj^^u 
end point of decay D2) can be seen more clearly over the distributions from Xa decays D4) 
and D5), which have higher mj-^^u edges. 

In Figs. 7c) and d), we nevertheless see some continuous background near the end point of 
mj-^^u which cannot be explained by Xa contamination. Note that for our choice of parameters, 
is considerably lighter than for the case studied in 
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moreover. 



is not too small 

compared to mq. The probability that one of the two hardest jets does not come from primary 
squark decays should therefore be higher than in the example analyzed in ref. 

' mj,u < 



These mis- reconstructed events also contaminate the m™^^ edge if we demand 



500 GeV < mj^u, as can be seen in Figs. 8a), b). Here we plot the higher of the two 
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Figure 7: a), b): mj^u distributions for a) mu < 55 GeV, b) 55 GeV < mii < 125 GeV. c), d): The same 
distributions after requiring rrij^u < 500 GeV < rrij^u. 
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Figure 8: Distribution of the higher of the two mj^i values for a) < mu < 55 GeV, b) 55 GeV< m/; <85 
GeV or 95 GeV < rrij-^i < 125 GeV with the cuts nij^u < 500 GeV< rrij^u. Figure c) corresponds to b) without 
the cuts on rnj^u and rrij^ii. 
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values in each event. However, the edges seem to be higher than the expected values in 
Table 3. Note that we exclude events with ma ~ niz because xt ~^ ^xt followed by Z — >• // 
has a higher rrij^i edge. For the sample with ma < 55 GeV, the contamination is seen as a 
change of slope, while for the samples with 55 GeV< ma < 85 GeV or 95 GeV< ma < 125 
GeV, no structure can be seen near the expected end point. 

This contamination actually was to be expected, because the events that fall above the real 
mj-^^a edge must be mis-reconstructed events where the jet originates from another sparticle 
decay or QCD radiation. Therefore the corresponding mji has no need to respect m^i^^, either; 
it tends to have a value larger than this nominal end point. The artificial upper limit of the 
mj-i^a distribution imposed by the cut then distorts the event distribution in Fig. 8 b). We 
find that the mj^i distribution without the requirement mj-^^a < 500 GeV < mj^a reproduces 
the m^^^ end point of decay chain D2) better for the events with 125 GeV> ma > 95 or 85 
GeV> ma > 55 GeV (Fig. 8c), although the distribution is still affected somewhat by events 
coming from decays. 

The lower edge m°^™ may be reconstructed from the mj^u distribution requiring mu > 
see Fig. 9. This distribution is much harder than the corresponding e/i distribution; 
this is a sign that the observed lower edge is real. The fit of the end point distributions will 
be given elsewhere . 




We now discuss the possibility to identify xt decays through the chain D2). In this case 
most daughter xt^s would decay further as Xi ^ ^ Xi; see Table 2, producing a r lepton in 
the last step of the cascade decay. The tU invariant mass never exceeds m -+ — m^o . Hadronic 
T decays might be identified by looking for a narrow jet that is isolated from other jet activity. 
Instead of studying the jet selection, we use information from the event generator to chose jets 
consistent with the parent r direction in the event list. The jet with minimum dR is selected 
as r-jet if dR < 0.3, \r]\ < 2.5 and Pr/PTj > 0.9, where P = Pr - Py^. We plot the mf^^ 
distribution, where jV is selected so that mj^a is minimal if the event contains several r— jets. 
When we compare the distributions for ma < 55 GeV (Fig. 10a) and 55 GeV < ma < 125 
GeV (Fig. 10b), we find the latter events clustered in the region mj^a < 190 GeV, while no 



17 



ID 23162 

UDFLVV 0.01)0 

OVFLVV 7.000 

ALLCHAN 0.4773E+05 




myiOGeV 



o 
o 

cu 
> 



ID 23163 

UDFLW 0.000 

OVFLW 3.000 

ALLCHAN 0.1275E+05 



™ myiOGeV 



a) b) 
Figure 10: The invariant mass distribution mj^u for a) mu < 55 GeV and b) for 55 GeV< mu < 125 GeV. 



such structure is found for the events with ma < 55 GeV. The only possible interpretation 
would be that most // pairs with 55 GeV< mu <125 GeV stem from the decay of a charged 
particle, xt- 

In the above plot, we are assuming 100% acceptance of r jets and no contamination from 
QCD jets. A rejection factor of O(IO^) against QCD jets together with a 40% r identification 
efficiency might be possible in the ATLAS experiment for jets with Pt > 30 GeV [Q. In 
[ p^ , the fake tau distribution is studied assuming a rejection factor of 15 for the case where 
X2 decays dominantly into r+r~Xi- Fake r backgrounds are then sizable in the region above 
the edge of the signal r+r~ distribution. The use of the jrll distribution to clean up the xt 
sample might nevertheless help to reconstruct edges from decay chain D2). 



4.2 Parameter Fitting 

In the previous subsection, we checked if it is possible to reconstruct the end points of invariant 
mass distributions involving charged leptons. Statistically, it seems possible to do so for decay 
modes Dl) and D2). This constrains mass differences among qi, X25 Xi? Cfl (Dl) and g^, xt , 
xt , vl (D2). We expect that these masses can be reconstructed with 0(10) GeV errors, as 
was the case in ref. [0. However, the corresponding errors on some MSSM parameters are 
significantly larger. 

In order to illustrate this point, we list two sets of MSSM parameters which reproduce all 
kinematic end points within Ax^ = 1, where Ax^ is defined as 



Ax' = ^(Mf - )V(AM,; 



(10) 



Here Mj runs over all five end points, m™f , 



m 



max(min) 



and m™*^^, of the decay chains Dl) 



and D2) listed in Table 3. We assume A(Mj) is 1% of M™^"* for distributions involving a jet, 
and 0.3 % of m™^^'™^"*. In Table 4, we list the solution with maximal and minimal yU (for 
tan/3 < 20) that satisfy Ax' < 1. 
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tan j3 






/X max 


575.381 


108.279 


196.758 


235.711 


20.0 


147.366 


202.577 


ji min 


554.638 


99.578 


218.696 


180.096 


20.0 


135.18 


185.826 



Table 4: Maximal and minimal /x solution satisfying Ax^ < 1 
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"^x2 
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11 max 


101.79 


163.15 


244.85 


284.94 


160.44 


285.10 


/X min 


89.0 


150.97 


190.41 


268.72 


144.19 


268.74 



Table 5: Neutralino and chargino masses in GeV for the maximal and minimal /j, solution. 

Note that the errors of the dimensionful parameters are strongly correlated, so that solu- 
tions with Ax^ < 1 almost fall onto a one-dimensional line in the seven-dimensional parameter 
space. Tabic 4 shows that the kinematic quantities wc have used in the fit give rather weak 
constraints for M2 and /j, with errors of order 20 GeV to 50 GeV. In fact, for fixed tan/? we 
find two distinct sets of solutions, with /x > M and /x < M, respectively. Moreover, one cannot 
fix the actual value of tan (3 from this fit; one can only determine that tan (3 ^ 8.65 where the 
minimum is achieved when M2 ~ 

Table 5 shows that the corresponding chargino and neutralino masses only vary within 
15 GeV between the two extreme solutions (except for Xg, which is almost not produced in 
q decays). Hence one will need additional information, beyond the kinematics of the decay 
chains Dl) and D2), to reduce the errors on the fundamental parameters. 

Reducing the errors on /j, and tan f3 would be necessary to predict the thermal relic density 
accurately. The fi (max,min) solutions predict ^x"^"^ ~ 0.160 and 0.122, respectively, as 
compared to 0.152 for the input point. The /x max point predicts a similar relic density as the 
input point; indeed, within the region with Ax^ < 1, we were not able to find solutions with 
fl^oh"^ > 0.165. On the other hand, the /i min solution predicts a significantly smaller relic 
density, and even smaller values of ^xoh"^ are possible if we relax the upper bound on tan/?, 
which was imposed "by hand" in this fit. For example, there is a solution with tan/? = 36 
and (Mi,M2,/x,mg,me«,mp) = (108.3,194.7,239.3,576.9,148.0,203.9) (all masses in GeV), 
giving fi^o/i^ ~ 0.112. We hence need to reduce the errors on both ji and tan/3. The former 
determines the size of the higgsino components of the LSP, which begins to be significant in 
this region of parameter space. The product yUtan/3 determines the amount of tl — tr mixing, 
which reduces the predicted relic density through a reduced fi mass and enhanced 5*— wave 
annihilation. In the following we discuss strategies that might be useful for reducing the errors 
on these two quantities. 
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Table 6: Squark branching ratios for the maximal and minimal (x solutions. 
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fi max 


fi input (univ) 


fi mill 


xl en 


2x 0.06 


2x 0.148 


2 X 0.105 


xt 


0.063 


0.096 


0.109 


Xl 


0.433 


0.402 


0.380 



Table 7: Branching ratios for the maximal and minimal /i solution. Wc assume universal soft sfcrmion masses 
and At = a,t the weak scale. 

One possibility is to measure some branching ratios. In Table 6, we compare the q decay 
branching ratios into charginos and neutralinos for the two solutions. Note that the ratio of 
the xt ^'^d X2 modes increases by more than a factor of three, from 0.45 to 1.71 for ul decay 
and from 1.11 to 4.25 for di decay, when switching from the /i max solution to the /i min 
solution. This is almost entirely due to the change of /i; the value of tan/? is not important 
here (as long as tan^/? ^ 1). 

The relative strengths of the signals from decay chains Dl) and D2) should thus yield 
important information to reduce the errors on MSSM parameters. The strengths of these 
signals can be extracted purely kinematically, e.g. from the relative number of events with 
mil below the X2 xt edge, respectively, and/or by trying to determine the fraction of 
di-lepton events that have a r— jet near the charged lepton pair, as discussed above. For a 
given solution in Table 4, all chargino and neutralino mixing angles and masses are fixed. As 
stated above, this is a fairly constrained fit where all relevant sparticle masses are effectively 
described by one parameter. The acceptances should then be very well calibrated from the 
mass constraints, so that systematic errors should be small. 

In order to extract squark branching ratios from the number of events with a lepton pair 
in the final state, one must know Br{xl —>■ cr) and Br{xt ~^ ^l)- These branching ratios 
also depend on MSSM parameters, as shown in Table 7. Here we assume that the f and 
i)r soft breaking mass parameters are the same as for first and second generation sleptons. 
We compute the f mixing angle by setting = at the weak scale; whenever it is sizable, 

— fji mixing is anyway dominated by the contribution oc /itan/5. With these assumptions 
all parameters required to compute these branching ratios can in principle be extracted from 
the kinematic fitting described above. 

The least critical quantity in Table 7 is the branching ratio for z/g xt decays. It decreases 
slightly with decreasing /z, due to the shrinking wino component of xt- However, this effect 
is weaker than the simultaneous increase of Br{xt ~^ ^e), which is due to the increasing wino 
component of xt- The strength of the signal from decay chain D2) is proportional to the 
product of these two branching ratios, which varies between 0.027 and 0.041. Together with 
the simultaneous change of Br{qL — > xt) shown in Table 6, this means that for our choice 
of parameters the signal strength of D2) decreases rapidly with increasing /i. Moreover, the 
relevant branching ratios do not depend significantly on the details of the f sector, and can 
thus be predicted fairly reliably from the quantities listed in Table 4.0 

Unfortunately this is not true for the branching ratio for x^ — > e/je decays, which does 
depend strongly on the mass and mixing angle of fi. Note that the prediction in Table 7 
for the input point (0.296) differs from the input value in Table 2 (0.236). This is because 
we ignored the reduction of soft breaking f masses through RG effects when computing the 

®In principle xt branching ratios could change somewhat if m^, < — to -+ , in which case additional 

2~body decay modes of xt would open up. However, such scenarios are already tightly constrained by LEP 
data, and would give rise to a variety of Higgs signals at the LHC. 
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entries of Table 7. In the given case these effects only reduce m^^ by ~ 5 GeV. This is sufficient 
to increase the partial width for X2 ~^ ^1'^ significantly; note that X2 ~^ decays are pure 
P-wave in the limit — 0, and the available phase space for these decays is not large in our 
case. Since for both fit solutions shown in Table 4 tan (3 is significantly larger than the input 
value, these solutions predict even lighter fi states and enhanced tl — tr mixing. The use of 
these parameters would therefore underestimate the true branching ratio for X2 ~^ ^R^ decays 
significantly. 

In order to extract the squark branching ratios of Table 6 to better than a factor of 2 
one will therefore need additional information on the f sector. This might be obtained by 
studying X2 ~^ 't^t'~x\ decays. As mentioned at the end of Sec. 3, it should be possible 
to determine the edge of the di— r— jet invariant mass distribution to ~ 5%. This would be 
sufficient to detect large differences between f and e masses; however, it would not suffice to 
distinguish between the three cases used in Table 7. To this end one would need to determine 
the ratio of branching ratios for x^ — >■ e^e~x\ and X2 ~^ t^^'X^ decays. The precision of this 
measurement might be limited by systematic effects, since the two signals have very different 
efficiencies.^ However, given that this ratio of branching ratios differs by more than a factor of 
4.5 between the three scenarios of Table 7 we think it likely that its measurement will help to 
reduce the errors of the extracted squark branching ratios significantly. Finally, once a linear 
collider of sufficient energy becomes available precision studies of fi production and decay will 



be possible |29, 30 



Let us summarize this somewhat complicated discussion by turning the argument around. 
One should first extract information about the f sector, e.g. by comparing signals from 
X2 — ^ T^^'X^ to those from X2 ~^ ^~^^~Xv This will give information on the soft breaking 
masses in the f sector as well as on the product jj, tan (3. This information, together with the 
result of the kinematic fit described above, will allow one to predict the branching ratios of 
the decays listed in Table 7 with reasonable precision. This in turn will allow to translate 
the measured strengths of the signals from decay chains Dl) and D2) into squark branching 
ratios. Finally, these branching ratios can be used to greatly reduce the error on fi. 

Another way to further constrain the relevant MSSM parameters is to include X4 decay 
edges from decays D3) and D4) in the fit. Just measuring values of these decay modes 

with 1% errors allows to reduce to 214 GeV. However, this would still not allow us to 
give an upper bound on tan (3.^ Moreover, this edge may not be visible for larger gluino and 
squark masses, where the production cross section is substantially smaller. 

Given that the very weak upper bound tan f3 < 20 which we imposed in the fit summarized 
in Table 4 is sufficient to predict Q^^^h"^ ~ 0.135 ± 0.03, it seems certain that the strategy 

outlined above will again allow to predict the thermal relic density to better than 20%. The 
only loophole occurs if tan f3 is very large. In this case the Xi ~ ^1 mass difference becomes 
so small that the r from fi decays become effectively invisible at hadron colliders. At the 
same time fi — Xi co-annihilation reduces the predicted LSP relic density by up to a factor 
of ten [^. One would then need to increase m and/or M in order to get a cosmologically 
interesting value of Q^oh'^; X2 ~^ ^R^ decays may not be open. In this case straightforward 
kinematic fitting as we described here will not be possible at the LHC, although one should 
still get a hint for the relative ordering of fi and cr masses by observing r— jets in missing 
i?T+jets events, which will yield the most robust SUSY signal in this case. In such a somewhat 
contrived scenario kinematical precision measurements would probably only be possible at a 



^Since the X2 — mass difference is quite small, one may have to allow one of the r— jets to be quite soft, 
"'^"including this new information gives the strong upper bound tan/3 < 11 at the la level, but tan/3 = 20 
remains allowed at 2a. 
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lepton collider. 



5 Discussion 

In this paper, we argued that LHC experiments can play a substantial role in predicting 
the contribution of thermal relic LSPs to the mass density of the Universe. Previous 
simulations in the literature were mostly done using mSUGRA assumptions, where usually 
l/ip ^ M^. In such a case the measured X2 cascade decay determines mg^, m^o and m^o. 
Once we know that the LSP is mostly bino, these three masses are sufficient to determine Q^^o 
within 20%. 

On the other hand, if Nature does not respect universality of all scalar soft breaking 
masses, it is possible that fi ~ M. In such a case, one needs to know fi and M in addition to 
m^o and mg^ to determine ^^^0, because s-channel exchange of Z and h, or WW production 
might play important roles in pair annihilation in the early Universe. We showed that, 
if yU ~ M, xt cind x1 will be produced copiously in qi decays. One can then determine all 
MSSM parameters needed to predict fl^A through the study of these cascade decay channels, 
using fits of kinematic end points and edges of invariant mass distributions. The isolation of 
xt ^ ^ Xi decays by observing the subsequent Xi ^ ^ decay may be used to improve the 
reconstruction of the xt production and decay kinematics. Moreover, the measurement of the 
relative number of events from different decay chains further constrains MSSM parameters. 
In the end it should again be possible to predict the thermal LSP relic density with an error of 
20% or better even in this more complicated scenario. In fact, this scenario is advantageous, 
since it allows us to determine both the gaugino and higgsino components of the LSP; these 
are needed to predict the strength of the LSP couplings to Higgs bosons, which in turn are 
required for predicting the LSP-nucleon scattering cross section. In mSUGRA scenarios with 
\fi\ > M one can probably only establish an upper bound on the higgsino component of the 
LSP, which only allows one to derive upper bounds on LSP-Higgs couplings. 

Notice that we only used information that can be extracted from studies at the LHC to ar- 
rive at this rather optimistic conclusion. If any one of the relevant masses could be determined 
with better precision elsewhere, e.g. at a lepton collider, the allowed region would shrink signif- 
icantly, since the fit of hypothetical LHC data resulted in an almost one-dimensional Ax^ < 1 
domain. 

In this paper we discussed the case where m <^ M so that neutralino decay into sfermion 
is open. This is a good assumption if Xi is gaugino like, as x? density overdose the Universe if 
m ^ M. For increased higgsino component of x?? such a requirement is no longer necessary. 
It is interesting to see if one can extract sparticle masses from decay distributions for such 
cases. 

In this paper, we did not study the case where fi <C Mi,M2, where the higgsino-like 
states Xi, X2 xt are nearly degenerate in mass. In such a case we should observe X2 > Xa 
production from squark decays in addition to Xs production, which is now mostly B. If scalar 
masses are not too large so that decays of the heavier neutralinos and charginos into real 
sfermions are open, the analysis is similar to the one that has been given in Sec. 4. If the 
sfermion decay mode is closed, the decay to a (virtual) Higgs boson might play an important 
role, unlike the case where /i ~ M. While on-shell Higgs bosons produced in SUSY cascade 
decays can be identified p4|, the kinematical fitting would be more difficult since it would be 



entirely based on jets. Note, however, that the thermal relic density of higgsino-like LSPs is 
small unless m^o ^ 500 GeV, in which case it might be difficult to even discover supersymmetry 
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at the LHC. 

We also did not discuss the case where M2 -C Mi,/i suggested in models with anomaly 
mediated supersymmetry breaking [Q. These models predict a rather heavy gluino. This 
results in a limited number of events even at the LHC, making precision studies rather difficult. 
In models with a not too heavy gluino while keeping M2 <^ Mi, /x, the relative number of events 
from B ^ cr and B might be useful to show that Mi 3> M2, if both modes are open. 

This would be sufficient to show that the thermal LSP relic density is small, independent of 
the relative ordering of M2 and /x, since both wino-like and higgsino-like LSPs with mass in 
the (few) hundred GeV range annihilate efficiently. 

We thus conclude that whenever LHC experiments find a large sample of SUSY events, it 
will be possible to either predict the thermal relic density of LSPs with a fairly small error, 
or else one will be able to conclude that thermal relic LSPs do not contribute significantly to 
the overall mass density of the Universe. In the latter case one would need physics beyond 
the MSSM, and/or a non-thermal LSP production mechanism, to explain the Dark Matter in 
the Universe. 
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